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Abstract. We consider the N-vortex problem on a ellipsoid of revo- 
lution. Applying standard techniques of classical perturbation theory 
we construct a sequence of conformal transformations from the ellipsoid 
into the complex plane. Using these transformations the equations of 
motion for the N-vortex problem on the ellipsoid are written as a formal 
series on the eccentricity of the ellipsoid's generating ellipse. First order 
equations are obtained explicitly. We show numerically that the trun- 
cated first order system for the three-vortices system on the symmetric 
ellipsoid is non-integrable. 
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1. Introduction 

The equations of point vortices motion in the plane were introduced by 
Helmholtz [1] and described as an Hamiltonian system by Kirchhoff in 1876 
[2]. These equations were first generalized to describe the motion of point- 
vortices on a sphere by Bogolmonov in 1977 [3]. Bogomolnov's work was 
put into solid mathematical grounds by Kimura and Okamoto in 1987 [5j. 
Independently Hally wrote in 1979 [4J the equations for vortices motion in 
symmetric surfaces of revolution (which englobe the sphere). To write his 
equations Hally represented the surface M on the complex plane with a met- 
ric conformal to the Euclidean metric: Let the i-th vortex, with vorticity Fj 
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be represented by the complex coordinate ^i. Calling the surface's conformal 
factor Hally's equations are 

(1) = h-\^n, &.) ( ^ -ij^ + ^te 

The difficulty in working with these equations is that one must know ex- 
plicitly the conformal factor , ^) and this factor is not known for many 
surfaces but the sphere. 

The aim of this paper is to write the equations for the N-vortex problem 
on the ellipsoid of revolution 

as an e-series when |e| << 1. To achieve this goal we write a perturbative 
coordinate transformation from the ellipsoid into the plane (see (|16p )and im- 
pose conformality at each order. It turns out that the transformation's n-th 
order term is given by the solution of a certain linear differential equation 
and that the first order term can be computed explicitly (see ([3])). Us- 
ing Hally's equations we calculated rigorously the first order perturbation 
term for the vortex equations on the ellipsoid. In spherical coordinates the 
truncated first order equations are given by the hamiltonian system with 
hamiltonian 

(2) H = Ho + €Hi 
where 

Ho = - ~ 2 cos(6'i) cos(6'j) - 2 sin(6'i) sin(%) cos((/)i - (f)j))] , 



1 ' 

Hi = - Y^r.Fj [(cos(0,)2 + cos(0,)2)] 

id 

and symplectic two-form 

sin{9i)d9i A 



N 

W 



The paper is organized as follows: In section ([2]) we introduce all the 
concepts used in the derivation of the vortex equations. Hally's equations 
are written as an invariant hamiltonian system and, for the spherical case, 
are shown to be equivalent to Bogolmonov's equations. In section ([3]) we 
derive the expression for the first order perturbation term. The main idea 
is to write the conformal transformation from the ellipsoid into the plane as 
a formal series in the eccentricity of the generating ellipse of the ellipsoid. 
In section ([4]) we exhibit two simple applications of equations ([2j). The case 
N = 2, despite being integrable, seems not to admit closed solutions. We 
show that, when \Ti\ = \T2\, the only solutions for which the two vortices 
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maintain their relative distance constant are the relative equilibria. For the 
A'^ = 3 case, we show numerically, by choosing specific parameters, that the 
three-vortices problem ([2]) on the ellipsoid is chaotic. 

2. The Equations of Hally 

In this section we discuss equations ([T]) from a geometric and invariant 
point of view. Hally's equations were deduced under the topological con- 
straint that the sum of all vorticities is zero and, in fact, this is essential 
in recasting the equations in Hamiltonian form. We also show, that, under 
this zero vorticity constraint, Hally's equations on the sphere, coincide with 
Bogolmonov's equations. 

2.1. Hally's equations as an hamiltonian system on the conformal 

plane. Let M be a 2-dimensional Riemannian manifold with metric tensor 
g. M is a symplectic manifold with symplectic form w given by the area 
form induced by g. Let : M — > M be the hamiltonian function. Hamilton's 
equations are 

(3) w{Xh ,-) = dH 

where Xh is the hamiltonian vector field induced by H. Now we assume that 
M is conformal to the plane, that is, that there are coordinates : M — > 
; i = 1, 2 where the metric is given by 

g {xi,X2) = h^{xi, X2) {dxi dxi + dx2 dx2) 

for some function h . In these coordinates the symplectic (area) form is 
given by 

w = h?'{xi , X2) dxi A dx2 . 
Introduce complex coordinates , ^) on the plane (xi , X2) by 

i= ^{xi+ix2) , 
i = ^{xi-ix2) . 

Using (C, we obtain 

g{i, = h\i, i) {d^ csd^ + ^ d^ 

and 

w = -ih^{C,C)dC AdC- 
Hamilton's equations ([3]) become 

h\^,Od^/\d^{XH,-)=d{iH) . 
For convenience we think of h^iC, d^ A dS, as the new area form and i H 
as the new hamiltonian function. Writing Xfj = i-§^ ~'~ equations 
become 

. o , - - dH 
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We now show that Hally's equations can be written as a hamiltonian 
system on the conformal plane. The equations are 

Multiplying each term of the sum in the right hand side by — 1^ and 
d 

adding the zero iVn log(/i(^fc , ^k)) to the right-hand side we obtain 



a = /i-'(?n,u(Efe'-^rfe^iog|c„-6P+ 



Assuming that the sum of the vorticities is equal to zero we have 

hHCn,QC = (E/rfcg|^iog|e„-efeP+ 

that is 

(6) h\^n,^n)^n = -i^! J-log {hi^n,^n) hi^k,^kMn - • 

k ^" 

But those are Hamilton's equations 

(7) w{Xh, •) = d{iH) 
for the symplectic two-form 

N 

(8) w = Y,^ih\^i,^i)d^iAd^i 

i=l 

and Hamiltonian function 

(9) H = ^ Fk r^log {h{^n, C"n) h{^k, ^kMn " 61') • 

k , n 

2.2. From Hally's to Bogolmonov's equations. Denote by S'^ C the 
two-dimensional sphere 

S^ = {ix,y,z)eW'\x^ + y'' + z'' = R^} . 

S'^ is a Ricmannian manifold with the metric induced by the Euclidean 
metric in M^. Introduce spherical coordinates on 5^ trough the relations 

X = Rcos{(/)) sin(0) , 
y = i?sin((/)) sin(^) , 
z = R cos{9) , 
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with 0<(^<2 7r, 0<e<7r. 

Consider the conformal map cr : — ^ C, a{9 ,(/))= u + iv where 

u = tan(^/2) cos(0) , 

(10) 

V = tan(6'/2) sin(^) . 
The conformal factor is given by 

where ^ = u + iv . Now let n , r2 € S'^ C M^. The Euclidean distance 
||ri — 'r2p in can be computed in the conformal variables. In fact 

||rL-r2f = {xi - + {yi - y2)'^ + {zi - Z2)'^ , 

= i?2 {2 - 2 cos(6li) cos(6'2) - 2 sin(6li) sin(6'2) cos(0i - (^2)} , 
= 4 cos(6li)2 cos{e2fR^ {(^^1 + + {ul + vD -2uivi-2u2 



i+€i€i 6 '^21 ) 



that is 



||ri - T2f = M6,6) ^(6,6) ICi - 61' ■ 

Now consider the map F : 5^ x • • • x — given by 
F{fi , r2 , • • • , rW) = {(T{ri) , a{r2) , cr(rAr)) 
The pull back by F of the symplectic form 

N 
i=l 

is just 



TV 



1=1 

where by Wi we denote the standard area form 

(11) Wi = Xi dyi A dzi + yi dzi A dxi + Zi dxi A dyi 

induced by the Euclidean metric on Sf. The pull back of the Hamiltonian 
function 

^=\Yl^k Tnlog {Kin. in) K^k, Un - Cfef ) 



2 

k , n 



IS 



(12) F*H = \Y^Tk r„iog {Wn - fuf ) . 



2 

k ,n 
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Hamilton's equations for Hamiltonian function ()12l) and symplectic form 
(jlip are precisely Bogolmonov's equations for N interacting vortices on a 
sphere. 

We notice that Rally's equations are equivalent to Bogolmonov's equa- 
tions only for the case where all the vorticities sums up to zero. Therefore, 
in dealing with the general vorticity case, one must use Bogomolnov's equa- 
tions. 



3. Equations on the Symmetric Ellipsoid - Perturbation 

Approach 

Let E"^ C represent the two dimensional ellipsoid of revolution 

Where R > and |e| << 1 . If e > the ellipsoid is prolate. In this case the 
eccentricity e of the generating ellipse is given by e = ^J^^^ ■ If e < the 

ellipsoid is oblate and the eccentricity of the generating ellipse is given by 
e = ■ We observe that e = O(e^) and therefore the truncated first order 
system is expected to approximate the real system for quite high values of 
the eccentricity. 

We introduce coordinates on E"^ trough the relations 

X = Rcos{(j)) sm{6) , 
(14) y = i?sin(0) sin(6') , 

z = i?^/r+7 cos(6l) , 

with 0<(/)<2 7r, O<0<7r.Ife = O the surface is a sphere of radius R. In 
this case the projection 

u = tan(6'/2) cos{(p) , 

(15) 

V = tan(6'/2) sin((/>) 

is a conformal transformation. We want to find a conformal transformation 
from E'^ to the conformal plane when e 7^ 0. By the symmetry of the 
ellipsoid the conformal factor must depend only on the radial variable r and 
not on the angular variable ^/J . We look for a transformation of the form 

= 2 arctan(r) + 2 E£,/i(r)eS 

(16) 

We observe that when e = the inverse of the above transformation is 
just transformation (llSp written in polar coordinates. 

The symmetry of the ellipsoid implies that the equator is a circle of radius 
R. In the {9, (j)) coordinates the equator has parametrization given hy 6 = ^ 
and (/> G [0, 2 vr) . Therefore in the conformal coordinates (u, v) the equator 
becomes the circle of radius r = + v'^ = 1 . We remark that the conformal 
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factor for the ellipsoid, when restricted to the circle of radius 1 is just the 
usual one for the sphere since the ellipsoid E"^ and the 5^ coincide at the 
equator. This is of utmost importance for what follows. In fact, since we 
will impose conformality at each order this trivial observation implies that 
= for alH. This is the condition that will guarantee the uniqueness 
of our series expansion. 

Remark 3.1. Alternatively (and equivalently) , the uniqueness of the ex- 
pansion can be obtained by imposing that the symmetry condition h{r) = 
^ h{r~^) has to be satisfied at all orders. 

The algorithmic determination of the fi{r) is obtained by imposing con- 
formality at each order and follows the standard procedures of classical per- 
turbation theory: Truncate the series at order n , find the conditions to be 
satisfied, and them assume those conditions to find the term of order n+1. 
To find we impose that 



generates a conformal transformation up to first order in e. 
The line element for the ellipsoid is 

ds"^ = sin2(^) #2 + {i + e ^{9)) d9^ . 

We want to write ds'^ in terms of r and (p. To this goal it suffices writing 
sin(^) and d9 in terms of r. But 

sin(^) = sin (2 arctan(r) -|- 2 /i(r) e) , 



9 = 2arctan(r) -|-2/i(r)e, 



(j) = tp. 



= sin(2 arctan(r)) -|- 2 cos(2 arctan(r)) /i(r) e + 0{e^) 



2 sin(arctan(r)) cos(arctan(r))-|- 



+2 (cos(arctan(r))2 — sin(arctan(r))2) fi{r) e -|- 0{e^) 



= 2^ + 2 




l+r- 



/i(r)€ + 0(e2). 



This gives 



+ 



8r(l -r^) 



/i(r)e + 0(e2). 



(l+r2)2 



(l + r2)2 



Also 




Therefore 
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Finally 



as - (1+^2)2 



1 + 



2(l-r^) 



/i(r)e +0(e2) 



+ 



l+((T^ + 2(l + r 



2^ 9/1 (r) 
9r 



} 



This will be conformal to the plane if it is a multiple of dip'^ + dr^. We 
observe that if e = then the transformation is conformal with conformal 
factor \J i^^_^^yi - For e 7^ the only way this transformation can be confor- 
mal for all values of e is if the coefficient of e in the first square bracket is 
equal to the coefficient of e in the second square bracket, that is if 

' dr r ' (l+r2)2 

This is a linear differential equation for Its general solution (for r > 0) 

is equal to 

{j^ +< 



(17) 



h{r) 



(l + r2 



■ r . 



Now by remark p.ip we must have /i(l) = what implies c = — ^. Finally 
we have 

l-r2 
/i(r) = —r. 



(1 +r 



2^2 



Transformation (jl6p becomes 

= 2 arctan(r) + r 6 + ^£2 (0 

(18) 

= V'. 

The conformal factor is, to first order given by, 

4i?2 



that is 

We write 

(19) 

where 



h{rf 
h{r) 



(l + r2)2 
2R 



(l-r2)2 



(l-r2)2 
1 + 4. 



ho{r) 



(l+r2) \ (l + r2)2 

h{r)=ho{r) + hi{r)e + 0{e^ 

(1-r^) 



2R 



2\2 



and hi{r) = 2R- 



(l + r2) '^"^ '""^'^ (l+r2)3 • 

For future reference we observe that in terms of the variables 9 and 4> we 
obtain 

ho{r) = 2Rcos^{e/2) and hi{r) = 4: R cos{9) cos'^{e/2). 
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Remark 3.2. The above expansion can not be carried explicitly further. In 
fact, the equation for the second order factor will be 

df2 lr(r2-l) I_i3r2+46r''-ri0_26r6 + 9r8 _ 
^ + 2 r2 + l ^ 4 (1 +r2)5 

The general solution for this equation is given by 



/2(r) 



where c is an arbitrary constant. The above integral can not be computed 
using elementary functions. 

3.1. N- Vortex Equations: Using the conformal factor just found we will 
write the N-vortex equations on the ellipsoid of revolution up to first order 
terms in e. We use the coordinates {9, (j)) introduced trough equations (fHI) . 
The symplectic two-form becomes 

w = VTTe d<f> A d{cos{e)) . 

Observing that 

Vl + e # A d(cos(6')) {Xh , • ) = ci</> A d{cos{9)) (Vl + e Xh , ■) , 
we define the vector field X = \/l + e Xh- X satisfies Hamilton's equations 
(20) d(l)Ad{cos{e))(x,-)=dH, 



and the flow (j)s{x) of the vector field X^ is a time reparametrization of 
the flow (/)t(x) induced by Xh- In fact t = vT+~e s. We therefore consider 
equation ()20p with if and /i given by ([9]) and (119^ respectively . H is given 
by 

^ = ^ Y^r, In [(/lo(ri) + e /ii(r,)) (/io(r-,) + e hi{rj) ) |e. - Cjf + 0{e^)] 

id 

that is 



+0(e2) 



The Hamiltonian function H is a function on the variables {O^cf))- More- 
over, the symplectic two- form (|20p is the standard symplectic form over the 
sphere. Therefore, we write 



H = \ HiJ^i [In (2-2 cosiOi) cos(0j) - 2 sin(0,) sin(0j) cos(</.i - </.j)) 

+e (cos(e,)' + cos(ej)2)] +0(e2). 
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We conclude that, to first order in e the dynamics of the N vortices on the 
elhpsoid is determined by the hamiltonian system with Hamiltonian function 

H = h H'ij^i^j [In (2 - 2 cos(6li) cosiOj) - 2 smiOi) smiOj) cos{(|)^ - (pj)) 

+e (cos(e,)^ + cos(0j)2)] 

and symplectic two-form 

N 

w = Ti d(j)i A d{cos{9i)) . 

i=l 

Alternatively, in cartesian coordinates the above Hamiltonian system be- 
comes 



(21) i/ = i ^ r, r, log \\n - r,f + 1 E ^^'^^^ + ' 

with symplectic two-form given by 

n 

(22) tt! = ^ Tj {xi dyi A dzi + yi dzi A dxi + Zi dxi A dyi) . 

i=l 

4. Applications: The N=2,3 cases. 

4.1. Two- Vortex problem. When N = 2 the hamiltonian system given 
by dill) and reduces to 



(23) 



n = T2j^^^ + eT2nxziz.^ 



f2 = r,^r^^ + eT,f2XZ2Z. 



Doing the following change of variables 

= Ti rl - r2 , 

(24) 

c = Ti r"! + r2 . 

we obtain 

'^'''^ fF +^ (Ai(c2,tt;2)c X z + A2{cz,w^)w x z) 



w 



(25) 

c = e {A2icz,Wz)cx z + Ai{cz,Wz)w X z) ; 

where 

"'"2 rir2 ' "'~2 TiF^ ' 

+ - ^2 + ^2 - u;^ 

Ai = 7=, 7=, — and = 7=, — . 

ATI ATI 

Equations (j25p are singular when ai c = , i.e. ri = r 2 . When e = 
equations (p5]) reduces to a linear system and can be easily solved. Those 
are the well known solutions for the two- vortex problem over the sphere: In 
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a system of coordinates such that its z-axis is in the same direction as the 
center of vorticity c the solutions are given in the ri, r2 variables by 

fi{t) = cos {9) , ji sin {9) , || c|| + - 

(26) 

^2(i) = it; (~/^ ^o*^ (^) ' ~^ si'^ (^) ' II ^1 

where 

9 = -it + (t), 

with (j) a constant and 

^ ^ 2 2,/ ~ II -||\2 ' 

a2 //^ + (a2 - ai II c ||) 

with 



// = if^. + = w — and ^ > 0. 

We observe that all solutions are periodic and the relative distance between 
the vortices is a first integral. 

When e 7^ the center of vorticity vector c is not preserved . In fact, 
only it's z component Cz is preserved. This first integral implies, by di- 
mensionality that system is integrable. Since ||ri|| = ||r2|| = 1 we obtain 
that 

||c|p + ||i(;||^ = fci and c-w = k2 

are constants of motion. Despite being integrable system (j25p seems not 
to admit explicit solutions for a general initial condition. Their relative 
equilibria can be easily characterized in the cases |ri| = |r2| = T . Then we 
obtain that 

||a1|2 + \\wf = T and c-w = 0. 

First we characterize the solutions such that ||ri(t) — r2(t)|| is constant along 
the motion. Since 

II ^ ^ l|2 II ^|2 2 I II ^||2 2 

IFi - ^^211 = ||c|| Q?! + ||u;|| 02 

and since for|ri| = |r2| either ai = or 02 = it follows that ||c|p and 
WwW^ are constants for the relative equilibria. Therefore 

(28) c- S= eA2C- {w X z) = 0. 

A2 = implies that Wz is a constant and therefore Wz = 0. From the equa- 
tions it follows that c • {w x z) = 0. This implies that the only way (|28|) can 
be zero is if c- {wx z) = 0. We arrive at the result that in a relative equilibria 
||c|p, ||?y|p and Wz are constants and that the vectors c , w and z are always 
in the same plane. The vectors c and w will precess around the z-axis with 
constant angular velocity and keeping their relative orientation. It follows 
from the definition that fi and r2 will also precess around the z-axis with 
constant angular velocity and keeping their relative orientation. Therefore 
the only solutions for which ||ri(t) — ^2(^)11 are constants are precisely the 
relative equilibria. 
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We compute the angular velocity of the relative equilibria in the cases 
where |ri| = |r2|. 

Case Fi = If Fi = r2 = T then ai = and 02 = . Assuming that 
Cx = A cos(r2t), Cy = A sm(0,t), = B cos(rit) and Wy = B sin(J7t) we 
obtain that the angular velocity is given by 



c 



2 



w\\^ 1" I 2r2 \ ' ^ I ■ 

We observe that the perturbation factor in the angular velocity is propor- 
tional to —75- and therefore the deviation from the corresponding spherical 

relative equilibria should be bigger for smaller vorticities. We also observe 
that when e = then, taking as before ||c|| = |c^|, the above expression 
reduces to (p7|) . 

Case Fi = — r2: If Fi = — r2 = T then ai = T^^ and 02 = . Proceeding 
as before we obtain 



^ = 77^ + — -T \ Cz+W, 



2 _ ^2 



Cl|2 I \^-' ' ^ j ■ 

The perturbation term is equal in both cases and again, as e — > we will 
have w — > 7 . 

4.2. The 3- Vortex problem. The symmetry breaking when e 7^ is ex- 
pected to destroy the integrability of the three-vortices spherical problem. 
In fact, when e = the center of vorticity vector c = Fi ri + r2 ^2 + Ts rs 
is conserved. The integrability follows by dimensionality. When e 7^ 0, only 
it's z-component = Fi zi + F2 Z2 + F3 2:3 is conserved. That the integra- 
bility is destroyed is shown numerically in this section. We do not pursue 
a complete numerical study of the problem. This will addressed in a future 
work. 

To calculate the Poincare sections for the three-vortex problem we will 
introduce an appropriate coordinate system. First we introduce a cylindrical 
system of coordinates (zi, 22, ^3, (Ai) 4>2^ (ps) with G [0, 2 7r) and Zi G (—1, 1) 
for i = 1,2. The transformation is defined trough 



X,; = a/1 - zf COs( 



yi = \ I- zf sin( 



Zi — Zi ■ 

Since 



3 

w 

i=l i=l 



Fj Xi dm A dzi + yi dzi A dxi + Zi dxi A dyj = ^ Fj d(j)i A dzi , 
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this is a symplectic transformation and the equations in cyhndrical coordi- 
nates are given by 

: I dH ^ I dH ^ 

0i = — and Zi = -— —— tor ? = 1, 2 . 

To reduce the symmetry we define the coordinates 

qi = Tizi- Ts 23 , 

gs = Ti zi + r2 Z2 +T3Z3, 

r,, — 4>1 \ 4)2 _ >) (t>3 
— 3 3 ^3 

P2 = ^-2f + f, 
P3 = ^ + f + f • 

It fohows that 

dqi A dpi + dq2 A dp2 + dq^i A dp^ = Fi dzi A + r2 dz2 A (i(/)2 + Fa dzs A (i</)3 . 

Therefore the new system of coordinates forms a canonical system. We 
observe that g3(t) = Cz and therefore q'i{t) = ^ = 0. Then ^33 is a cychc 
variable and the Hamiltonian ()2ip is given by 

= H{qi,q2,pi,p2,Cz). 

Since this Hamiltonian is now defined in a four dimensional phase-space, 
the Poincare maps can be visualized. The equations are integrated numer- 
ically using the Fehlberg-7(8) Runge-Kuta. The relative error in energy is 
kept smaller then lO"^'^. A constant step size equals to 10^'^ is used in all 
integrations. 

Figure 1 show the pictures of three Poincare maps. The onset of chaotic 
behavior is clear. It appears (as expected) around the homoclinic solution 
present in the integrable case e = 0.00. We observe that for e = 0.05 
and e = 0.10 the eccentricity of the generating ellipses are 0.22 and 0.33 
respectively. 
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Figure 1: Poincare maps for the Three-vortices problem. Here Fi = r2 = = 
3, Energy = — |. Cj; = — |;. (^2 is the horizontal axis and p2 is the vertical axis. 
The values of e are, from top to bottom, e = 0.00 , 0.05 , 0.10 respectively. 
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